Simulation techniques for cosmological simulations

Modern cosmological observations allow us to study in great detail the evolution and history of the large scale structure hierarchy. The fundamental problem of accurate constraints on the cosmological parameters, within a given cosmological model, requires precise modelling of the observed structure. In this paper we briefly review the current most effective techniques of large scale structure simulations, emphasising both their advantages and shortcomings. Starting with basics of the direct N-body simulations appropriate to modelling cold dark matter evolution, we then discuss the direct-sum technique GRAPE, particle-mesh (PM) and hybrid methods, combining the PM and the tree algorithms. Simulations of baryonic matter in the Universe often use hydrodynamic codes based on both particle methods that discretise mass, and grid-based methods. We briefly describe Eulerian grid methods, and also some variants of Lagrangian smoothed particle hydrodynamics (SPH) methods.Comment: 42 pages, 16 figures, accepted for publication in Space Science Reviews, special issue "Clusters of galaxies: beyond the thermal view", Editor J.S. Kaastra, Chapter 12; work done by an international team at the International Space Science Institute (ISSI), Bern, organised by J.S. Kaastra, A.M. Bykov, S. Schindler & J.A.M. Bleeke

Magnetic Field Amplification in Galaxy Clusters and its Simulation

We review the present theoretical and numerical understanding of magnetic field amplification in cosmic large-scale structure, on length scales of galaxy clusters and beyond. Structure formation drives compression and turbulence, which amplify tiny magnetic seed fields to the microGauss values that are observed in the intracluster medium. This process is intimately connected to the properties of turbulence and the microphysics of the intra-cluster medium. Additional roles are played by merger induced shocks that sweep through the intra-cluster medium and motions induced by sloshing cool cores. The accurate simulation of magnetic field amplification in clusters still poses a serious challenge for simulations of cosmological structure formation. We review the current literature on cosmological simulations that include magnetic fields and outline theoretical as well as numerical challenges.Comment: 60 pages, 19 Figure

Particle acceleration mechanisms

We review the possible mechanisms for production of non-thermal electrons which are responsible for non-thermal radiation in clusters of galaxies. Our primary focus is on non-thermal Bremsstrahlung and inverse Compton scattering, that produce hard X-ray emission. We briefly review acceleration mechanisms and point out that in most astrophysical situations, and in particular for the intracluster medium, shocks, turbulence and plasma waves play a crucial role. We consider two scenarios for production of non-thermal radiation. The first is hard X-ray emission due to non-thermal Bremsstrahlung by nonrelativistic particles. Non-thermal tails are produced by accelerating electrons from the background plasma with an initial Maxwellian distribution. However, these tails are accompanied by significant heating and they are present for a short time of <10^6 yr, which is also the time that the tail will be thermalised. Such non-thermal tails, even if possible, can only explain the hard X-ray but not the radio emission which needs GeV or higher energy electrons. For these and for production of hard X-rays by the inverse Compton model, we need the second scenario where there is injection and subsequent acceleration of relativistic electrons. It is shown that a steady state situation, for example arising from secondary electrons produced from cosmic ray proton scattering by background protons, will most likely lead to flatter than required electron spectra or it requires a short escape time of the electrons from the cluster. An episodic injection of relativistic electrons, presumably from galaxies or AGN, and/or episodic generation of turbulence and shocks by mergers can result in an electron spectrum consistent with observations but for only a short period of less than one billion years.Comment: 22 pages, 5 figures, accepted for publication in Space Science Reviews, special issue "Clusters of galaxies: beyond the thermal view", Editor J.S. Kaastra, Chapter 11; work done by an international team at the International Space Science Institute (ISSI), Bern, organised by J.S. Kaastra, A.M. Bykov, S. Schindler & J.A.M. Bleeke

Search for gravitational waves from Scorpius X-1 in the second Advanced LIGO observing run with an improved hidden Markov model

We present results from a semicoherent search for continuous gravitational waves from the low-mass x-ray binary Scorpius X-1, using a hidden Markov model (HMM) to track spin wandering. This search improves on previous HMM-based searches of LIGO data by using an improved frequency domain matched filter, the J-statistic, and by analyzing data from Advanced LIGO's second observing run. In the frequency range searched, from 60 to 650 Hz, we find no evidence of gravitational radiation. At 194.6 Hz, the most sensitive search frequency, we report an upper limit on gravitational wave strain (at 95% confidence) of h095%=3.47×10-25 when marginalizing over source inclination angle. This is the most sensitive search for Scorpius X-1, to date, that is specifically designed to be robust in the presence of spin wandering. © 2019 American Physical Society

Search for Tensor, Vector, and Scalar Polarizations in the Stochastic Gravitational-Wave Background

The detection of gravitational waves with Advanced LIGO and Advanced Virgo has enabled novel tests of general relativity, including direct study of the polarization of gravitational waves. While general relativity allows for only two tensor gravitational-wave polarizations, general metric theories can additionally predict two vector and two scalar polarizations. The polarization of gravitational waves is encoded in the spectral shape of the stochastic gravitational-wave background, formed by the superposition of cosmological and individually unresolved astrophysical sources. Using data recorded by Advanced LIGO during its first observing run, we search for a stochastic background of generically polarized gravitational waves. We find no evidence for a background of any polarization, and place the first direct bounds on the contributions of vector and scalar polarizations to the stochastic background. Under log-uniform priors for the energy in each polarization, we limit the energy densities of tensor, vector, and scalar modes at 95% credibility to Ω0T<5.58×10-8, Ω0V<6.35×10-8, and Ω0S<1.08×10-7 at a reference frequency f0=25 Hz. © 2018 American Physical Society

Erratum: "A Gravitational-wave Measurement of the Hubble Constant Following the Second Observing Run of Advanced LIGO and Virgo" (2021, ApJ, 909, 218)

[no abstract available

MADRE: The Memory-Aware Data Redistribution Engine

We report on the development of a new computational framework for efficiently carrying out parallel data redistribution in a limited memory environment. This new library, MADRE (The Memory-Aware Data Redistribution Engine), is an open-source, C/MPI-based toolkit de-signed for quick and easy integration into application codes that have demanding data migration needs. At the same time, MADRE exposes a lower-level API that greatly facilitates the devel-opment and incorporation of new algorithms into the MADRE framework, thus serving as a potential organizing entity for continued research in this area. Finally, we develop, describe, and test in detail several new parallel redistribution algorithms that are incorporated into the MADRE distribution.